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ABSTRACT 

SN 201 lja was a bright (I = -18.3) Type II supernova occurring in the nearby edge on spi¬ 
ral galaxy NGC 4945. Flat-topped and multi-peaked Ha and H (3 spectral emission lines ap¬ 
pear between 64 - 84 days post-explosion, indicating interaction with a disc-like circumstellar 
medium inclined 30-45 degrees from edge-on. After day 84 an increase in the H- and K-band 
flux along with heavy attenuation of the red wing of the emission lines are strong indications 
of early dust formation, likely located in the cool dense shell created between the forward 
shock of the SN ejecta and the reverse shock created as the ejecta plows into the existing 
CSM. Radiative transfer modeling reveals both « 1.5 x 10 1 Mq of pre-existing dust lo¬ 
cated ~ 10 16 7 cm away and « 5 x 10 -5 Mq of newly formed dust. Spectral observations 
after 1.5 years reveal the possibility that the fading SN is located within a young (3-6 Myr) 
massive stellar cluster, which when combined with tentative 56 Ni mass estimates of 0.2 Mq 
may indicate a massive (> 25 Mq) progenitor for SN 201 lja. 

Key words: circumstellar matter — stars: winds, outflows — supernovae: general — super¬ 
novae: individual (SN 2001 lja) 


1 INTRODUCTION 

Type IIP, the most common type of core collapse supernovae (CC- 
SNe), have broad (~10 4 km s _1 ) hydrogen emission lines along 
with a near constant “plateau” of optical luminosity throughout the 
first ~100 days. The widely accepted progenitors of Type IIP SNe, 
red supergiants (RSGs), have masses ranging between q£- 25 
and mass loss rates of ~10~ 6 to 10 -4 Mq yr _1 dChevalier et al.l 
120061 : Mauron & Josselinffjolll) . Type Iln SNe also show broad hy¬ 
drogen emission, but, in addition, they show narrow (~ 100 km 
s _1 ) hydrogen emission due to ionization of the surrounding, dense 
pre-existing circumstellar material (CSM). These SNe likely have 
more massive progenitors such as Luminous Blue Variables (LBV), 
or Wolf-Rayet (WR) stars that can have mass loss rates which are 
orders of magnitude larger, 1 0~ 5 to 10 -2 Mq yr _1 dKiewe et al.l 


ers, anywhere from hours to years after explosion dLeonard et 


al. 


l2012l:ISmith et alEoi lLhoij) . Whether Type IIP or Type Iln, these 
massive star progenitors can undergo periods of dramatic mass loss 
prior to explosion which has a direct impact on the SN evolu¬ 
tion. Over the past decade we have seen numerous observational 
signatures of the SN ejecta interacting with previously shed lay- 
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2000l: Andrews et al.l201(jl : lFransson et al.l20 1 4l : I Smith et al.l2015l : 

Maeda et al.ll2015l . for example). By observing the SN-CSM inter¬ 

action we can get a better understanding of pre-supernova mass 
loss, pathways of dust production, and the link between progenitor 
and SN type. 

Due to their short lifetimes and ability to return material back 
to the ISM quickly, CCSNe are the likely cu lprits for the dust 
produ ction in dusty high-z galaxies dGall et alJl201 ll : iDwek etHll 
120141) . This is a double-edged sword, as they are also efficient 
destro yers of du st, at least in prese nt day galaxies iTemim et al.1 
120151) . Although lSlavin et alJ <1 20151) propose effective SN dust de¬ 
struction, but only within a subset of the appropriate assumed pa¬ 
rameters. Many recent studies of nearby CCSNe have searched for 
signatures of dust formation and estimated the dust masses. The 
results of these studies in dicate small amounts o f newly-formed 
dust, 10 -2 - 10~ 4 Mq 1 Elmhamdi et alJ [20031 : Sugerman et al.1 
120061: iMeikle et al.ll2007l : iKotak et al.ll2009l : lAndrews et alJl20ld. 
for example), much less than the 0.1-1 Mq needed to account for 
the excess of dust seen in the early Universe. Within the past few 
years, far-IR studies have revealed a few sources which may hold 
promise to unlocking the dust mystery. Herschel observations of 
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Figure 1 . Color composite Gemini/GMOS image of SN 201 lja (indicated in white). Optical standards (Table 3) are indicated by blue solid circles, NIR 
standards (Table 4) by dashed yellow circles. 


three nearby young supernova remnants indicate the presence of 
0.1 Mq of cool dust in Cas A d Barlow et alJl2Q10h . 0.4 - 0. 7 Mq 
of cool dust in the ejecta of SN 1987A l Matsuura et alj20 111,120151) 
and 0.18-0.27 Mm of cool d ust in the Crab Nebula I Gomez et al.l 
l2Q12l : IOwen & Barlowll2015l) . Recent SOFIA observations of a su¬ 
pernova remnant at the center of the Milky Way revealed 0.02 Mq 
of warm (100 K) dust, which also has seemed to have survived the 
passing of the reverse shock dLau et al~1l2015t ). This strongly sug- 
gests that almost all of the dust in SNe is form ed after > 1000 days 
^Wesson et alJl2QI5l : iBevan & Barlow! l2015h . While reservoirs of 
cold dust may help solve this problem, there is a growing body of 
evidence that grain growth may be occurring in the ISM enriched 
by these early SNe and could be the main site for dust production 
at high-z dMichalowskfco 1 5l ). However, a plausible mechanism for 
growing refractory grains in the ISM has yet to be identified. 

Observational signatures of dust formation in CCSNe man¬ 
ifest themselves in several different forms. The optical luminos¬ 
ity will decrease while almost simultaneously the NIR will in¬ 
crease as the dust grains absorb the shorter wavelength light and 
re-emit it in the IR. The grain formation will also alter the op¬ 
tical spectra, creating asymmetric and blue-shifted lines as the 
dust grains attenuate the red, receding side of the ejecta prefer¬ 
entially. This observational evidence of dust formation has now 
been seen in numerous SNe, including the nearby an d well-studied 
SN 1987A Jlucy et al.l 1 1 9891: I Wooden et al.l 1 1 9931) . SN 2003gd 
dSugerman et alJl2006l: i Meikle et alJl2007l) . SN 2004et dSahu et aid 
120061 : Kotak et al,l2009l) and many others. While it was initially be¬ 
lieved that the dust grains could only condense 300-600 days after 
explosion, when the ejecta had expanded and cooled, there have 
been more and more confirmed cases of dust forming much ear¬ 


lier, within ~100 days of explosion. This can occur due to shock 
interaction with nearby CSM creating an area between the forward 
and reverse shocks with temperatures and densities appropriate for 
grain growth, this area is known as the cool dense shell (CDS). For 
example SN 1998S showed dust formation signatures between days 
140 - 268 dLeonard et al.ll2000l) . and SN 2005ip appears to have 
formed dust both in the C DS between day 75-150 and then again 
in the ejecta after day 750 dSmith et all2009l : lFox et all2009l) . The 
bright Iln SN 2010j 1 sho ws c ontinuous dust formation between 40 
and 240 days dGall et al.l2014l) . Although not classified as Type Iln, 
the Type Ib/c SN 2006jc also form ed dust via CSM in t eraction be¬ 
tween 50 - 75 days post-explosion dMattila et al.ll2008l : ISmith et~ahl 
2008) and the Type IIP SN 2007od formed d ust sometime betwee n 
day 120 - 230 through the same mechanism dAndrews et al.ll2oToh . 
As we will present below, with more and more long-term moni¬ 
toring of CCSNe, there seems to be evidence of non-Type-IIn SNe 
exhibiting signs of CSM interaction, even months after explosion. 
This not only allows a separate channel for dust formation in CC¬ 
SNe, but can also reveal important properties of SN evolution. 

SN 201 lja was discovered in NGC 4945, an edge-on spi¬ 
ral (Figure 1) loca ted at a distance of 3.36 ± 0.09 Mpc 
dMouhcine et alj2005l) . NGC 4945 is one of the closest Milky Way 
analogs with a near solar metallically, particularly at increasing 
distances from the galactic center where the supernova is located 
dstanghellini et ali2015h . This is of importance, particularly in dust 
formation, as higher me tallicity galaxies tend to be able to pro - 
duce higher dust masses dDraine et all2007l : lGalametz et al,l201 ll) . 
On 2011 December 19, optical spectra obtained of the object indi¬ 
cated that it was a Type II supernova which most closely matched 
SN 2004et about a week after maximum dMonard et alfeOlU) . al- 
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though maximum optical light of SN 201 lja does not appear to 
have been achieved until ~ 20 12 Janu ary 1 4 (Fig ure 2). Radio ob¬ 
servations presented in IChakraborti et ahl ( 20131) . further suggest 
an explosion date of 2011 December 12 (JD 2455907). Therefore, 
throughout this paper we choose 2011 December 12 as day 0 for 
SN 201 lja, and day 34 as the date of maximum light. In Section 2, 
we discuss the data reductions and present a comprehensive analy¬ 
sis of these data in Section 3. Section 4 contains a discussion of the 
implications of the data, and Section 5 includes analysis in regards 
to dust formation, CSM location, and progenitor characteristics. Fi¬ 
nally, in Section 6 we briefly summarize the significant results of 
the paper. 


2 OBSERVATIONS AND DATA REDUCTION 

Optical imaging and spectra were obtained with GMOS/Gemini 
South (GS-2012A-Q-79, GS-2013-Q-93, PI Andrews). A Color 
composite Gemini/GMOS image of SN 201 lja is shown in Figure 
1, and a summary of observations and resultant photometry in Ta¬ 
ble 1. The gVi' images were reduced and stacked using the IRAfO 
gemini package. The instrumental g'r'i' magnitudes were trans¬ 
formed to standard Johnson-Cousins VRI using tertiary standards 
created from stars in the field (Fi gure 1 and Table 2) and transfor¬ 
mations presented in lWelch et al.l 12007 ). Uncertainties were calcu¬ 
lated by adding in quadrature the transformation uncertainty quoted 
in lWelch et all (' 2007 1. photon statistics, and the zero point devia¬ 
tion of the standard stars for each epoch. Figures 2 and 3 show the 
resultant light curves. 

For each GMOS/Gemini South epoch, three spectra of 900s 
were obtained in semester 2012A and six spectra of 900s were ob¬ 
tained in semester 2013A. The spectra were obtained in longslit 
mode using grating B600 and a slit width of (//75. With a resolv¬ 
ing power R = 1688, this corresponds to a velocity resolution of 
~ 180 km s -1 . Central wavelengths of 5950, 5970, and 5990 A 
were chosen to prevent important spectral features from falling on 
chip gaps. A 2x2 binning in the low gain setting was used. Spectra 
were reduced using the IRAF gemini package. The sky subtraction 
regions were determined by visual inspection to prevent contamina¬ 
tion from material not associated with the SN, and the spectra were 
extracted using 15 rows centered on the SN. The spectra from each 
individual night were averaged and have been corrected for the ra¬ 
dial velocity o f NGC 4945 (460 km s -1 ) obtained from narrow line 
emission dMonard et al 1201 Ih . Spectrophotometric standards were 
not taken as part of our Gemini program but are taken at least once 
per semester. For each of our observations we have flux calibrated 
our spectra using the spectrophotometric standard that is closest in 
date to the observation. The final spectra are presented in Figures 4 
and 5. 

Seven epochs of JHK S imaging, spanning 8 to 272 days 
post-explosion, were obtained from NTT/SOFI data made avail¬ 
able on the ESO archive (184.D-1140(N), PI Benetti). A summary 
of the epochs and resultant photometry is presented in Table 3 and 
shown in Figure 2. Observations were taken with a pixel scale of 
0.288"/pixel in a four-point dither pattern for background subtrac¬ 
tion. Images were corrected for crosstalk and background contam¬ 
ination using IRAF scripts provided on the SOFI website then re- 

1 IRAF is distributed by the National Optical Astronomical Observatory, 
which is operated by the Association of Universities for Research in Astron¬ 
omy, Inc., under cooperative agreement with the National Science Founda¬ 
tion. 


Table 1. Gemini/GMOS Photometry of SN 201 lja 


Day 

JD 

V 

R 

I 

84 

2455991 

17.20 ± 0.08 

14.87 ± 0.06 

13.76 ± 0.06 

112 

2456019 

18.00 ±0.10 

15.71 ±0.06 

14.62 ± 0.06 

159 

2456066 

18.37 ±0.05 

16.10 ±0.05 

15.03 ± 0.04 

450 

2456357 

19.88 ±0.10 

18.42 ±0.11 

17.64 ± 0.10 

508 

2456415 

19.95 ±0.10 

18.76 ±0.12 

17.88 ±0.10 

807 

2456714 

19.74 ±0.11 

18.61 ±0.11 

18.05 ±0.08 


Table 2. Tertiary VRI Standards for NGC 4945 


Star 

a (J2000) 

m h + 

3 (J2000) 
-49° ± 

V 

R 

I 

A 

.341633 

.495212 

16.30 

15.42 

14.98 

B 

.350917 

.506569 

15.39 

14.90 

14.73 

C 

.353678 

.507214 

15.70 

14.43 

14.93 

D 

.345111 

.517944 

16.31 

16.26 

15.67 

E 

.240586 

.501742 

16.44 

15.82 

15.57 

F 

.243756 

.514639 

16.60 

16.11 

15.56 

G 

.297883 

.493965 

15.85 

14.56 

15.12 

H 

.279450 

.516431 

16.28 

14.88 

15.45 

I 

.297790 

.494104 

16.77 

15.81 

16.02 


Photometry is from NOMAD and DENIS catalogs. 


duced, aligned, and stacked using standard IRAF techniques. Dif¬ 
ferential aperture photometry was performed using the 2MASS 
standards listed in Table [4] On 2011 December 17 (Day 15) and 
2012 April 30 (Day 113) we also obtained JHK S imaging from 
the 0.8m Infrared Imaging System (IRIS) telescope at Cerro Arma- 
zones Observatory, Chile (operated by the Astronomical Institute 
of the Ruhr-Universitat Bochum and the Universidad Catolica del 
Norte.) Aperture photometry was also performed using the stan¬ 
dards presented in Table 4. Additional low resolution optical spec¬ 
tra for days 64 and 237, observed with NTT/EFOSC, were also 
obtained from the ESO archive (184.D-1151(Z), PI Benetti). Ob¬ 
servations were taken using a slit width of l'.'O and grism 13, which 
has a resolution of ~ 2.8 A/pixel. Standard IRAF reduction proce¬ 
dures were applied, and background subtraction was done using the 
same area as for the GMOS spectra. Flux calibration was achieved 
using spectrophotometric standards taken on the night of observa¬ 
tion. 

Spitzer IRAC (3.6 and 4.5 /im) images were obtained at four 
epochs from day 105-857. The images were mosaicked and resam¬ 
pled using standard MOPEX procedures to improve photometric 
quality. A summary of the observations and fluxes are shown in 
Table 3. PSF photometry was performed using the Point Response 
Functions (PRF) developed for IRAC. Additional mid-IR obser¬ 
vations at 18.72 and 10.77 /im were obtained with VLT/VISIR 
(288.D-5031(A), PI Wesson) on 2012 February 26 (day 76) and 
March 14 (day 93), respectively . For the B10.7 filter a total 3600s 
was taken on target as was 7200s for the Q2 filter. Both observa¬ 
tions resulted in a non-detection, so only upper limits can be in¬ 
ferred. These are 2 mJy for 10.7 fim and 10 mJy for 18.72 /.im. 






















4 Andrews et al. 


Table 3. Near- and Mid- IR Photometry of SN 201 lja 


Day 

JD 

J 

H 

K 

3.6 fim (mJy) 

4.5 fim (mJy) 

8 

2455915 

11.46 ±0.05 

11.85 ± 0.03 

10.76 ± 0.06 



15 

2455922 

11.12 ±0.04 

10.93 ± 0.04 

10.32 ± 0.04 



40 

2455947 

10.78 ± 0.03 

10.29 ± 0.03 

9.92 ± 0.06 



65 

2455972 

11.15 ±0.03 

10.63 ± 0.02 

10.20 ± 0.07 



93 

2456000 

12.59 ± 0.03 

12.02 ± 0.04 

11.60 ±0.06 



105 

2456012 

- 

- 

- 

12.54 ± 0.29 

14.35 ± 0.24 

113 

2456020 

12.85 ±0.10 

12.22 ± 0.10 

11.60 ±0.15 



121 

2456028 

13.12 ±0.02 

12.45 ± 0.04 

11.87 ±0.06 



243 

2456150 

14.11 ±0.11 

12.65 ± 0.07 

11.49 ±0.07 



272 

2456179 

14.39 ± 0.04 

12.91 ± 0.02 

11.63 ±0.04 



486 

2456393 

- 

- 

- 

13.48 ± 0.29 

15.32 ± 0.24 

637 

2456544 

- 

- 

- 

10.16 ±0.26 

11.63 ±0.22 

857 

2456764 

- 

- 

- 

7.79 ± 0.20 

8.17 ±0.17 


Table 4. Tertiary 2MASS Standards for NGC 4945 


Star 

a (J2000) 
I96 h + 

d (J2000) 
-49° ± 

J 

H 

K 

1 

.265558 

.507378 

12.68 

12.44 

12.38 

2 

.279346 

.516434 

14.77 

14.35 

14.27 

3 

.264956 

.535126 

15.08 

14.64 

14.55 

4 

.316055 

.519047 

15.79 

15.01 

14.87 

5 

.287417 

.534851 

15.50 

14.90 

14.84 


3 ANALYSIS 

3.1 Internal Extinction 

The initial discovery spectrum, shown in Figure 4, suggested 
a larg e amount of extinction toward SN 201 lja IMonard et al.l 
120111). The Milky Way f oreground is only Ay =0.48 mag 
i Schlaflv & Finkbeinetl[2oTTI) so much of the extinction is inter¬ 
nal to NGC 4945, which is not unexpected as_ extinction estimates 
along its galactic plane are A„ > 11-13 dMarconi et al Jl2Q()ol) . In 
order to determine an accurate extinction value for SN 201 lja, 
we have employed a variety of tests. Because of the low reces¬ 
sion velocity of the galaxy (450 km s -1 ) the use of Na ID doublet 
was insufficient for this galaxy due to the blending of internal and 
external Na absorption lines. This is not necessarily a hindrance, 
as it is still under debate whether this method is acceptable for 
SNe ( Poznans ki et alJl20lH . l2012h . especially in the presence of 
circumstellar material. Therefore, we compared our unreddened, 
early-time spectra with optical spectra of SN 2004et, a prototypical 
Type IIP, from similar epochs and applied extinction corrections 
unt il the spectra were coincident. SN 2004et has an E(B-V) = 0.43 
dSahu et al.l2006l) . and comparisons on both day 7 and day 84 yield 
a total E(B-V)=1.8 in SN 201 lja, using the reddening law of CCM 
dCardelli et al J198(d) . Throughout this paper we use a value of total 
reddening, foreground to SN201 lja, to be E(B-V) =1.8. 

As a second method and sanity check, we measured the rela¬ 
tive line strengths of the Fla and FI (3 nebular emission, assuming 
both Case A and Case B recombination in the HII region surround¬ 
ing the SN, as these are the two extremes of hydrogen recombina¬ 
tion, and the area surrounding SN 201 lja must lie somewhere in 
between. In the absence of reddening we would expect a ratio of 
2.85 for Case B and 1.91 for Case A, but we find a ratio closer to 
17.25 suggesting 1.64 < E(B-V) < 2.01 dSeatordl 1 9791) . This is in 


reasonable agreement with comparison to the SN 2004et spectra. 
When we measure this ratio on day 159, E(B-V) has increased to 
2.10, indicating an additional increase in A„ of 0.55 mag between 
these two dates, assuming R„ = 3.1. As we will discuss below, this 
is likely due to early-dust formation. 

3.2 Optical Lightcurve Evolution 

Figure 2 shows the VRI evolution of SN 201 lja between days 28 - 
1138. The I-band AAVSO observations up to day 77 indicate that 
the SN may have been increasing in luminosity since discovery and 
reached a maximum brightness on 2012 January 14, day 34. Shown 
in detail in Figure 3, the plateau phase only lasted for ~ 35 days, 
with a decline rate of 0.01 mag day -1 . The drop into the radioactive 
decay portion of the light curve that begins on day ~70, lasted only 
15 days with a Al of 1 magnitude, which was similar for the other 
optical bands. It is possible that the transitional phase lasted until 
day 112, and if this is the case, the SN faded by 1.8 mags in 42 
days (a rate of 0.04 mag day -1 ). We have reason to believe, as 
explained below, that the start of the radioactive decay phase was 
in fact around the time of our first GMOS observation of SN 201 lja 
on day 84. 

Comparison with the lightcurves of other Type II SNe shown 
in Figure 3, particularly after day 50, shows that SN 201 lja may lie 
somewhere between a Type IIP and a Type IIL classification-wise. 
The blurring between “plateau” and “linear” designations of Type 
II SNe is becoming more common as the sample of well studied 
objects increases, with a survey of hydrogen rich SNe presented 
bv lAnderson et alj 120151) concluding that there are not two distinct 
classes, but rather a continuum based most strongly on the enve- 
lope m ass at explosion. This is somewhat at odds with lArcavi et al .1 
d2012l) who make a contrary claim that there are distinct groups 
based on progenitor types. The short plateau duration, the absolute 
magnitude at maximum (I = -18.3), and the steep drop into the ra¬ 
dioactive decay phase all point to a CCSN with a smaller hydrogen 
envelope, which we will discuss in more detail below. 

Assuming the radioactive decay phase started on day 84, the 
V-band magnitudes should then be 17.47 on day 112 and 17.92 on 
day 159 if the only energy source is the decay of 56 Co (dashed line 
in Figure 2). This is roughly 0.5 magnitudes brighter than the ob¬ 
served magnitude, requiring an increase in optical depth (r) of 0.5 
between day 84 and 112 to account for the magnitude deficit. Mod¬ 
eling of the optical and IR photometry lend credence to this hypoth¬ 
esis, and as we will show in Section 4, radiative transfer modeling 
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Days after Explosion 


Figure 2. Optical and IR light curves of SN 201 lja from values listed in Tables 1 and 3. The NIR curves have been shifted up one magnitude for clarity. V and 
I photometry prior to day 84 has been obtained from AAVSO. The dashed line follows the decay of 5 ®Co. 


shows that there is a substantial increase in tv between day 105 
and day 486. By itself this does not necessarily confirm an increase 
of extinction prior to day 159, but in tandem with the optical spec¬ 
troscopy presented in the next section, we have sufficient reason to 
believe new dust was being formed within 100 days of explosion in 
SN 201 lja. 

Alternatively if the radioactive decay phase is assumed to be¬ 
gin on day 112, the SN does not suffer from any increased extinc¬ 
tion and behaves as expected until ~400 days, when the luminosity 
levels out through our last observation on day 857. This late-time 
plateau can be caused by various physical mechanisms, notably a 
scattered-light echo, shock-interaction, or ambient light from the 
parent stellar cluster. These scenarios are discussed in detail in Sec¬ 
tion 5.2. 


3.3 Spectra 


and red-shifted peaks are seen in Ha at -1400 km s , and 2900 
km s^ 1 , in addition to the narrow nebular line at 0 km s _1 . There 


H /3 emission shows peaks at -2100 km s _1 and 3100 km s -1 , but 
not the existence of a second red-shifted component. We cannot as¬ 
certain if a double-peaked structure exists in the He I A5876 line 
due to the blending with the Na ID doublets from local and host 
galaxy extinction. Multipeaked asymmetric hydrogen lines have 
been seen in other SNe, such as SN 1993J dMatheson et al.| |200C ), 
SN 1998S (Leonard et alj|2000h. S N 2004di dChugai et alj|2007l>. 
SN 2007od ^Andrews et al.l 1 201 ol> . SN 2009ip jMauerhan et al.l 
|2014|) . and PTF1 liqb dSmith et al.l2015l t. and are mostly attributed 
to a torrodal or disc geometry of surroundin g CS M material. Com¬ 
parison with models presented in iFransson et al.l d2005h indicate 
that the CSM around SN 201 lja may be represented by a torus 
inclined between 30-60 degrees with an angular thickness between 
30 and 50 degrees. 


Early-time spectra of SN 201 lja show a prototypical Type II su¬ 
pernova (Figure 4, Day 7), with blue-shifted Ha and H/7 emission 
lines, peaking at -3500 km s -1 and -3000 km s -1 , respectively, 
and very little else. While most other SNe show a very blue con¬ 
tinuum at these early times, the high reddening towards SN 201 lja 
represses this feature. This is further exemplified by a strong Na ID 
absorption feature present around 5900A. A high-velocity feature 
(HV) may also be present in the Day 7 spectrum at ~13,700 km 
s -1 , likely due to flash-heate d unshoc ked pre-explosion mass loss 
from the progenitor star dChugai et al .1120071) . 

Our first GMOS spectra of SN 201 lja was obtained on day 
84, when the multi-component hydrogen lines that seem to appear 
in the day 64 spectra are now prevalent (Figure 5). Blue-shifted 


Other broad nebular lines such as Ba II A6142, [Sc II] A5527, 
A5658, and A6246, and [O I] AA6300.6363 are dominant in the 
nebular spectra. Fe II AA 4924,5169 may also be present and some¬ 
what blended with H/i and other Sc lines. An intermediate-width 
component of [O III] A 5007 A seems to appear at this time as 
well, and persists until our last spectra on day 807. Alternatively, 
this could be a -735 km s -1 blue-shifted peak of Fe A 5018 A. If 
this is the case, then the red emission peak of H j3 is likely blended 
with a similarly blue-shifted component of Fe A4924 A. The only 
other published example _of broad [O III] emission is in the well- 
studied Type Iln 1995N Lransson et alj|2002h . In that object, lines 
of a similar width (1500 km s _1 ) were seen, the only difference be¬ 
ing in SN 1995N they showed a much more boxy profile, similar to 
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Figure 3. Optical I-band lightcurves for a sample of Type II SNe. The ri ght figure sh ows a zoomed-in regi on of the firs t 200 days. Data are from: SN 
1987A jHamuv et alJll988l : Tsuntzeff et~al]| 19881: lHamuv & Suntzerfll990h . SN 2007it jAndrews et afll201ld) . SN 2004et JSahu et alj|2003) . SN 1999em 
jElmhamdi et all2003 ), and SN 2002hh jPozzo et al .12 006) 



Figure 4. Optical spectra of SN 201 lja. Day 7, 64, and 237 are from NTT/EFOSC2, and the rest are from Gemini/GMOS. All spectra have been flux calibrated 
but shifted by a constant for presentation. An extinction correction of E(B-V) = 1.8 has been applied to all spectra. 


what we see in our Her emission lines. Those authors attributed this 
emission to unshocked ejecta composed of oxygen core material 
coupled with a low hydrogen envelope mass. 

By days 159 and 237, when the SN is well into the nebular 
phase, the strong P-Cygni profile in Ha has all but disappeared. 
In addition to having multiple peaks, the Ha (and to some ex¬ 
tent H/3) lines begin to show signs of flattening as the nebular 
phase persists. Boxy, flat-topp ed spectr a have also been see n in SNe 
1993J dMatheson et alfcOOOli . 19 98S jLeonard et al .11200(1) . 2004et 
dKotak et al.ll2009l) . and 2007od dAndrews et alJl20loh . and are at¬ 
tributable to CSM interaction. Also noticeable in the Ha emis¬ 


sion is the degradation of the red peak from the day 64 and 84 
observations (Figure 5). Over the course of ~180 days, the ratio 
of the blue and red peak emission has evolved from 1:1 to 2:1. 
Further comparison between Ha emission of SN 201 lja at ~day 
84, SN 2004dj at day 128, and SN 2007od at day 232 (Figure 6) 
shows similar structures, particularly in the peaks at ±1500 km 
s -1 . This double-peaked structure was attributed to the interaction 
of the ejecta with a torus or disc of material surrounding the su¬ 
pernova, with the stronger attenuation on the red side due to dust 
forming in the line of sight. The extra peak at +3200 km s^ 1 is 
more of a mystery. It is possible that the multiple peaks seen in the 





































SN201 lja Evolution 7 




Figure 5. Ha; (top) and H/3 (bottom) evolution for a subset of epochs from 
Figure 4. The enhanced red emission in H (3 on day 84 could be due a blend 
with Fe A 4924 A. 

hydrogen lines could come from asymmetries in 56 Ni in the ejecta. 
This non-uniformity could cause uneven ionization and excitation 
in the ej ecta, as was suggested for the He I lines in the Type lib SN 
2008ax jTauhenberger et ahll201 lh . 

The Ha line width of 5000 km s -1 persists from day 64 un¬ 
til our last observation on day 807. Late-time spectra show the 
emergence of strong lines of [N II] AA 6548,6583 and [S II] AA 
6716,6731 Afsee Figure 4). More striking is the blue continuum 
that emerged in the last two epochs. This occured at the same time 
as the leveling off of the optical lightcurve discussed in section 3.2. 
The two are likely caused by the same phenomenon, and are dis¬ 
cussed in detail in Section 5.3. 


3.4 IR Observations 

In Figure 2 we also show the first ~ 300 days of JHK evolution, 
including a maximum around day 40. Between day 8 and 40 there 
is an increase in the absolute magnitude of all NIR bands, although 
the largest increase of 1.5 mag occurred in the H-band. In the weeks 
following the maximum there is a steep decline, very similar to the 
optical lightcurve evolution, until day 121. Between day 121 and 
243 the K-band brightens by ~ 0.4 magnitudes and an almost con¬ 



Figure 6, Ha profiles of oth er SNe with k nown CSM i nterac tion. Data 
are from lAndrews et aD j2010h (2007od) and Ivinko ct al.l J2006I) (2004dj). 
All three SNe were classified as normal Type II, yet all three show the 
same double-peaked emission profile, with a dominant blue-peak within 
4-8 months post-explosion. 


stant plateau in the H-band. This corresponds to the same time pe¬ 
riod when we believe there was an increase in optical fading by 
0.5 magnitudes, and that the Ha emission showed the most red- 
side attenuation. A similar increase in the H- and K-bands was 
observed around da y 50 in SN 2005ip dFox et ai] 12009 ) and SN 
2006jc dSmith et al.ll2008h . In both cases, the increase in the NIR 
can be attributed to dust grain formation. Although there is no NIR 
photometry available to us after day 272, the brightness is once 
again declining into a radioactive tail, indicating that the new grains 
have either been destroyed in the reverse shock or have cooled suf¬ 
ficiently to only be detectable at longer wavelengths. 

3.5 56 Ni Mass Estimates 

We estimate a 56 Ni mass for SN 201 lja of Mjvi= 0.22±0.03 Mg 
employing the methods of lHamu a®, using the V magnitude 
at day 84. The bolometric luminosity of the radioactive tail (in erg 
s -1 ) is calculated as 

, . — [Vt — A(V) + BC] + 5logwD — 8.14 

log w ^t = --, ( 1 ) 

2.0 

with Av =5.6 mag, and a bolometric correction of BC = 0.26. The 
nickel mass is then calculated at various times during this phase, as 

tt —tp _ 

M m = {7M6xlQ- AA )L t exp[ *+* 26 ' ]M 0 . (2) 

Here 6.1 days is the half-life of 56 Ni and 11.26 days is the e- 
folding time of 56 Co and tt-to is the age of the SN. This calculated 
Ni mass is much larger than other Type IIP SNe. There is, of course, 
some uncertainty in the extinction correction, not to mention the lu¬ 
minosity at each epoch could be underestimated due to additional 
dust formation or overestimated due to circumstellar interaction. 
If this is an accurate measurement of the 56 Ni mass, it would put 
SN 201 lja in a regime of high mass normally reserved for stripped 
envelope CCSNe, particularly lib and lb ( ILvman et al.ll20141) . For 
comparison, other “normal” Type Up SNe such as SNe 1999em, 
2003gd, and 2004dj each have 56 Ni masses ~0.02 M 0 , or a full or- 
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der of magnitude lower than estimated here I El mhamdi et al.l2003l : 
iHendrv et alJl2005l:IVink6 et all2006l) . 


4 RADIATIVE TRANSFER MODELING 


To determine the amount and type of dust being formed in SN 
201 lja, we have used the MOCCASIN 3D Monte Carlo radia¬ 
tive transfer code fercolanqjh_aL[ |2005}._and references therein). 
As was done for SN2007it l Andrews et al.ll201 lah and SN 2010jl 
dAndrews et al.ll201 lbl ) we assume three different dust geometries. 
The “smooth” distribution describes dust distributed uniformly 
within a spherical shell surrounding the SN. The “torus” model dis¬ 
tributes dust uniformly within a torus at some inclination around the 
SN. Finally, the “clumpy” model places dust in clumps and scatters 
it throughout a spherical shell surrounding the SN. For grain sizes 
we are using a standard MRN grain size distribution of a~ 3 ' 5 be¬ 
tween 0.005 and 0.05 /xm ( iMathis et all 19771 ). For each epoch, we 
use the Spitzer IRAC data to set the date, and extrapolate the op¬ 
tical and NIR photometry from the surrounding observations. We 
also correct all photometric points for the assumed E(B-V)=1.8. 

Continuing on in the manner presented in Andrews et al. 
(2010, 2011a,b), the dust and luminosity for the source was lo¬ 
cated between an inner radius Ri„ and an outer radius R otl t of 
a spherically expanding shell, with the diffuse emission luminos¬ 
ity being proportional to the density at each location. The smooth 
model assumes the density of dust in the shell was inversely pro¬ 
portional to the square of the radius. For the clumpy model the pho¬ 
tons originate in the inhomogeneous interclump medium, where the 
clumps are considered to be optically thick and spherical. For the 
torus models, densities are specified for the inner and outer walls, 
with the dust distribution falling off linearly between the two radii 
fercolano et al .1 120071) . For each model we used grains of amor - 
phous carbon (AC) using the optical constants of Iflannerl dl988t) . 
due to the lack of rnid-IR detection at 10.8 /xm, which would 
have indicated a strong silicate component. Inputs for the models 
were luminosity, ejecta temperature, inner and outer radii, and dust 
masses. For each epoch and parameter, initial estimates were ac¬ 
complished using blackbody fits to the optical and IR data (solid 
curves Figure 7). This yielded dust temperatures (Td) of roughly 
550 K for the first epoch (day 105), and 725K for the remaining 
epochs. These temperatures are consistent with typical warm dust 
temperatures. Inputs and outputs for each epoch and geometry are 
listed in Tables 5 and 6. 

For the shell models, the first epoch (day 105) is best fit by an 
Ri„ of 5.3 x 10 16 cm, which roughly corresponds to the evapora¬ 
tion radius of the initial flash of the SN. Epochs 2 and 3 (day 486 
and 637 respectively) used an Ri„ of 4.5 x 10 15 cm, and for our 
final epoch (day 857) Ri„ = 3.8 x 10 15 cm. For the torus models, 
the best fits were achieved with slightly different values for Ri„. 
On day 105 we assumed Ri„ = 5.0 x 10 15 cm, and for the remain¬ 
ing epochs R;„ = 9.0 x 10 15 cm. We also kept the temperature of 
the ejecta (T e /) at 7000K for all epochs, since this is a reasonable 
ejecta temperature l lSugerman et al.l2006h and does an adequate job 
fitting the first epoch of optical points. The combination of strong 
Ha and [O I] emission lines in the R band relative to the continuum, 
and at later times the additional blue flux has made it unlikely that 
the visible photometry is consistent with any reasonable blackbody 
temperature. Keeping those two values constant, we then varied the 
luminosity, dust masses, and for the torus models the number den¬ 
sity along the inner edge to get the most accurate fits. For all epochs 
we found that the clumpy distribution predicted on average about a 


factor of two higher dust masses than the smooth and torus models. 
We have also found that varying the grain sizes and exponent on the 
grain distribution alters the dust mass outputs. In general increas¬ 
ing the minimum grain size does not alter the dust mass, increasing 
the exponent integer on the grain size distribution will decrease the 
dust mass, and increasing the maximum grain size from 0.05 /tm 
to 1 /xm will increase dust mass by roughly a factor of two. These 
values are within the range of masses created by varying the dust 
geometries (Table 6), which is our largest source of uncer tainty. 
A recent in-depth study presented in iBevan & Barlow! J2015h finds 
that a MRN distribution could not reproduce the line profiles of SN 
1987A at early epochs, but as the grain size distribution is not our 
largest source of uncertainty we retain our standard MRN distribu¬ 
tion for consistency. 

We estimate a dust mass for our first epoch of 1.0 x 10 -4 Mq 
for the smooth model, 1.9 x 10 ~ 4 A/q for clumpy, and 1.4 x 
10~ 4 Mq for the torus. The remaining epochs yielded smooth dust 
masses of 6.7 x 10~ 5 Mq for day 486, 6.0 x 10 -5 Mq for day 637, 
and 3.8 x 10~ 5 Mq for day 857. The larger dust mass from the first 
epoch is likely due to the flash heating of pre-exisiting CSM dust 
that is just outside of the shock radius, and is roughly 200K cooler 
than the newly formed dust in following epochs. On day 105, the 
maximum distance the shock could have traveled assuming a veloc¬ 
ity of 11000 km s' 1 is 1 x 10 16 cm making this scenario entirely 
plausible. Additionally R i„ for this first epoch is larger than R ou t 
for the remaining epochs, indicating a different location of the mod¬ 
eled dust. This dust will cool and the IR echo will fade over the next 
year, so that by the time of Epoch 2, the SED should be dominated 
by newly formed ejecta dust. 

Around the time of Epoch 2, the ejecta should have traveled 
an average distance of 2.5 x 10 16 cm (assuming a late-time expan¬ 
sion of 6000 km s _1 measured from Ha widths.) This places the IR 
emission, and therefore the newly formed dust, mostly in the CDS 
interior to the point of CSM interaction. Of course this is just the SN 
envelope that is traveling at this velocity. The metal-rich core mate¬ 
rial itself is likely moving much slower, closer to the 1500 km s _1 
seen in SN 1987A l lMcCravil993ll . and will only be at a distance of 
6.3 x 10 15 cm. The steady decline of estimated dust masses over 
all of the epochs is most likely due to the dust in the distant CSM 
responsible for the IR echo cooling, leaving only the newly formed 
dust as the main near-IR contributor. It is also possible that as the 
reverse shock propagates in Epochs 3 and 4 it destroys the dust that 
was formed in the CDS between Epochs 1 and 2, or that some of 
the newly formed dust has cooled as well. This is fully consistent 
with the picture presented in the spectral evolution. 

On day 105, the MOCASSIN models predict ar„ = 0.03. Us¬ 
ing methods presented in IFox et all 12009 ). we can also indepen¬ 
dently estimate the optical depth. If we assume L max = 1.25 x 10 9 
Lq for SN 201 lja based on the lightcurve maximum, then E ~ 
1.25 x 10 49 erg using the estimates of SN 2005ip and our previ¬ 
ous estimation of SN 2007it. We can constrain the IR energy using 
the black body fits, which results in an Lbb of 5.6 x 10 6 Lq and 
therefore an IR energy over the first 105 days of 1.97 x 10 47 erg. 
Inserting these numbers into r ~ E ^ EE , we calculate ar = 0.015, 
consistent with our MOCASSIN fits. 

An increase in t v between epochs 1 and 2, as mentioned 
above, occurs coincidentally with the emission line and optical 
lightcurve extinction, and therefore the formation of new dust. 
While the total dust mass is actually smaller in the second epoch, 
the volume into which it is contained is much smaller as well, since 
R OU ( has decreased by over a magnitude between the two, requir¬ 
ing similar amounts of dust in a much smaller volume. This region 
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Figure 7. SEDs of days 105 (upper left), 486 (upper right), 637 (lower left), and 857 (lower right). All photometric points have been corrected for E(B-V)=1.8, 
and optical and NIR fluxes have been extrapolated from other observations to match the Spitzer observation dates. Blackbody fits are shown as solid lines. 
Mocassin fits are shown in dashed (smooth), dotted (torus, 45°), and dash-dotted (clumpy). Note the y-scale is different for Epoch 1 than for the other Epochs. 


Table 5. Monte Carlo Radiative Transfer Shell Models 


Epoch 

T ej (K) 

Rin (cm) 

Rout (cm) 

Ltot. (L 0 ) 

Tv 

Smooth 

M d (M 0 ) 

Clumpy 

M d (M 0 ) 

105 d 

7000 

5.3el6 

5.3el7 

1.8e8 

0.03 

1.0e-4 

1.9e-4 

486 d 

7000 

4.5el5 

4.5el6 

1.3e7 

2.60 

6.8e-5 

1.4e-4 

637 d 

7000 

4.5el5 

4.5el6 

1.0e7 

2.34 

6.1e-5 

1.2e-4 

857 d 

7000 

3.8el5 

4.2el6 

9.1e6 

1.76 

3.9e-5 

8.6e-5 


Table 6. Monte Carlo Radiative Transfer Torus Models 


Epoch 

T ej (K) 

Rin (cm) 

Rout(cm) 

Ltot. (L 0 ) 

Face-On (0°) 

Tv M d (M 0 ) 

T v 

45° 

M d (M 0 ) 

Edge-On (90°) 

Tv M d (M 0 ) 

105 d 

7000 

5.0el6 

5.0el7 

1.8e8 

0.0 

1.3e-4 

0.03 

1.5e-4 

0.05 

1.8e-4 

486 d 

7000 

9.0el5 

4.0el6 

1.3e7 

0.0 

4.8e-5 

0.0 

6.1e-5 

3.31 

7.7e-5 

637 d 

7000 

9.0el5 

4.0el6 

9.8e6 

0.0 

4.3e-5 

0.0 

5.4e-5 

2.94 

6.8e-5 

857 d 

7000 

9.0el5 

4.0el6 

8.9e6 

0.0 

3.5e-5 

0.0 

4.3e-5 

2.37 

5.5e-5 
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is also almost entirely interior to the original R e „ ap , again pointing 
to newly formed grains. In Section 3.2 we estimate an increase in 
A„ of 0.5 between day 84 and 112, model results would therefore 
indicate that the most likely dust geometry would be in a torus in¬ 
clined at roughly 45°. As can be seen in Table 6, between face-on 
and 45°, r = 0, but increases from 0.5 to 3.31 for inclinations be¬ 
tween 47-90°. Overall, the total mass of new, warm dust formed 
in SN 201 lja appears to be ~ 6.0 x 10~ 5 Mq by day 857, which 
although consistent with dust masses of other CCSNe, is still con¬ 
siderably small er than the amount needed to account for the dust 
seen at high-z dMorgan & Edmundsll2003b . Although in the Intro¬ 
duction, we have described the much larger dust masses seen at 
longer wavelengths for SN1987A and other SNe observed at much 
later epochs. 


5 DISCUSSION 

5.1 Circumstellar Interaction and Environment 

The expansion of the supernova shock into the surrounding envi- 
ronment can create X-r ay emission if the CSM is dense enough. 
Ichakraborti et al.l 120131) found that the X-ray flux in SN 201 lja in¬ 
creased by over a factor of 4 between day 29 and day 113 which was 
attributed to the supernova interacting with a CSM created from a 
non-steady wind. As discussed in detail above, the optical spectral 
evolution implies that CSM interaction is occurring by day 84 and 
per sists onward, agreeing w ith the X-ray observations. In particu¬ 
lar, [Qiakraborti_et_alj (|201J) hypothesize the progenitor underwent 
a blue supergiant (BSG) phase shortly before explosion, producing 
a low density cavity immediately surrounding the SN. Other CC¬ 
SNe showing an increase in X-ray emission over the first 120 days 
include the type IIP SN 1 999em (Pooley et al. 2002) and the type lb 
SN 2006jc (Immler et_al.||2008[), and both instances were attributed 
to shock interaction with the CSM. 

If we use the optical spectra to constrain the date of the onset 
of shock interaction, we can safely assume that while it was already 
beginning on day 64, it was fully occurring by day 84. Assuming 
the initial expansion velocity of 11 000 km s -1 was sustained, the 
CSM is located 8 x 10 15 cm away from the SN. On day 113, when 
the X-ray observation was taken, the shock had reached a maximum 
distance of 10 16 cm. 

It is important here to remark on the evaporation radius 
(Rejjap), the cavity cleared by the explosion of the supernova. Fol¬ 
lowing the evaporation radii computed bv lFransson et al J J2014I ) for 
SN 2010jl, we have estimated the ranges of R e „ap from grain size 
and type. Using an absolute magnitude of I = -18.3, and standard 
bolometric correction of 0.26, we estimate L ma x = 4.8 x 10 42 erg 
s _1 , or l/6 th the maximum luminosity of SN 2010jl. This means 
that for graphite grains (which are used in our MOCASSIN mod¬ 
els), and a T e // = 10 4 , R e „ ap = 2.4 x 10 16 and 9.0 x 10 16 cm 
for 1.0 fim and 0.001 fim, respectively. Therefore we can reason¬ 
ably assume an R e „ ap = 5 x 10 16 cm. If of course the BC or 
E(B-V) estimations are under- or over- estimated this radius can 
change. A lower L ma x of 1 x 10 42 erg s _1 could bring R et) ap as 
close as 10 16 cm , which was th e R e ^qp estimated for SN 2005ip 
Jfqx et al J|20 loi) and SN 2007it dAndrews et alj|2011al ). Addition¬ 
ally, a clumpier CSM could also allow grains to survive closer to 
the SN. Therefore, R e „ ap is consistent with the CSM interaction 
beginning by day 84. 

As mentioned above, an intermediate-width [O III] emission 
line of ~ 3000 km s -1 may be visible in our day 84 spectra and per- 
sist up to the last spectral observation. According to lFransson et al .1 


(20021) . if the explosion is asymmetric it is possible to get [O III] 
A5007 emission at early times. The scenario described for a simi¬ 
lar feature seen in the Type II 11 SN 1995N is that the intermediate 
component is actually the oxygen core of the supernova. The pres¬ 
ence of this line also may help explain the interesting shape of the 
hydrogen lines. To get high enough ionization energies an asym¬ 
metric explosion, with most of the energy occurring in the plane, 
is needed. This would create a slower expansion velocity along the 
plane and much faster in the polar direction. This scenario could ex¬ 
plain the intermediate-width Ha and H/3 lines as well, which have 
an underlying broad component coming from the polar region. 


5.2 New Dust Formation and Uocation 


Over the past few years it has become apparent that there are two 
modes of dust formation in CCSNe. If there is strong CSM inter¬ 
action occurring it is possible for dust grains to condense in the 
cool dense shell (CDS) between the forward and reverse shocks. 
This has been seen in other SNe such as SN 1998S (Pozzo et al 


2004) , SN 2004dj jMeikle et al J201 lb. SN 2005ip dFox et al.l200' 


s 

a 


Smith et al.l l2009b. S N 2006ic dSmith et al.1 l2008h, SN 2007od 
([Andrews et all201Ch . and SN 2010jl dSmith et all2012l : IOall et all 
12014b . Taken solely on the formation of grains within the first 150 
days of explosion in SN 201 lja would point to dust formation in 
the CDS. The appearance of asymmetrical hydrogen lines in the 
optical spectra is a smoking-gun. We cannot rule out the possibility 
of small amounts of dust forming in the expanding ejecta at later 
times, but observational evidence points to new dust forming be¬ 
tween day 84 - 112 in a region between 4.5 x 10 15 and 4.5 x 10 16 
cm away. 


5.3 Fate-time luminosity 

After day 400 there is little change to the luminosity of SN 201 lja 
in VRI. There are three possibilities to this late-time brightness: a) a 
scattered light echo off the surrounding CSM. b) radiative shocks, 
or c) fading of the SN into the ambient brightness of the parent 
cluster. These three possibilities are not mutually exclusive. Scat¬ 
tered light echoes occur when light from the supernova scatters off 
circumstellar or interstellar dust and redirects it back into the line 
of sight of the observer at some later time. The scattering is more 
efficient at shorter wavelengths, causing the light echo spectrum to 
look bluer than the peak SN light. At the same time, a leveling out 
of the optical lightcurve is created as the original light pulse trav¬ 
els outward through the intervening material creating a constant 
luminosity. The late time spectra, SED, and lightcurve all show in¬ 
dications that this could be a possibility, as the spectrum becomes 
bluer and the light curve flattens out. The inference of pre-existing 
dust from the SED modeling as well as the high level of extinction 
in the direction of SN 201 lja also imply that there is plenty of in¬ 
tervening circumstellar and interstellar dust to create light echoes. 
Unfortunately, when one creates a light-echo spectrum of the SN 
using the integrated fluency from the spectra prior to day 100, sim¬ 
ilar to that done for SN 1980K in lSugerman et al.l d2012h while the 
emission lines appear similar, there is not enough flux in the blue to 
create a match (Figure 8). So while a light echo may be present, it 
would be much fainter than the late-time spectrum presented here, 
so we have ruled this out as the cause of the late-time luminosity. 
A caveat to this assumption is that the SN explosion date could be 
earlier than estimated here, which would mean the inclusion of an 
earlier, bluer light echo spectrum. Deep HST imaging of the region 
around SN 201 lja could resolve a scattered light echo if one exists. 









































SN2011ja Evolution 


11 



Figure 8. Day 807 spectrum (purple) of SN 201 lja showing enhanced blue emission. Comparison with the light echo spectrum created from an integrated 
fluency of the first 84 days (red) indicates that a light echo cannot be responsible for the flux blueward of 6000A . The orange and yellow spectra are synthesized 
stellar populations created with Starburst99 for 3 and 6 Myr. It is possible that the late-time luminosity has a large component of the parental stellar cluster. 


The second scenario requires dense CSM to convert the ki¬ 
netic energy of the expanding supernova into radiation via shocks. 
This is a fairly common occurrence in Type Iln SNe at early times, 
and seems to be more common than originally observed in nor¬ 
mal Type II as observations are extending past the first 3 months 
after explosion. By day 807, the bulk of the SN ejecta will have 
reached a distance of 4.1 x 10 16 cm (or 7.7 x 10 16 for the fastest 
moving ejecta). This places the edge of the ejecta clearly outside 
of Reuap, and at the distance of the CSM responsible for the early 
IR luminosity, making radiative shocks plausible. Unfortunately a 
shock driven luminosity, while theoretically being responsible for 
the constant light curve, would not make the optical spectrum bluer. 
Additionally, the CSM should slow the ejecta expansion, and there 
does not seem to be a substantial decrease in the expansion speeds 
in Ha from day 450 to day 807. Therefore we assume that radia¬ 
tive shocks are not responsible for the late time optical plateau, 
although additional X-ray observations could potentially reveal if 
CSM interaction is indeed occurring. 

Since only ~5-10% of massive stars are formed in isolation 
dOev et alj|2013l) , it would not be out of the question to detect a 
CCSN still located in the stellar cluster from which it is formed that 
would become detectable as the SN faded. For instance, the type IIP 
SN 2004dj occurred in the massive cluster Sandage-96 and showed 
very similar late-time photometric and spectroscopic evolution as 
SN 201 lja dVinko et al.ll2009l ). The lightcurves leveled out after 
day 800 to the magnitude of the cluster in pre-explosion HST im¬ 
ages. The spectrum also showed a very blue continuum below 6000 
A, but the nebular features of the SN (Ha, [O I]) are still present in 
the red part of the spectrum. Inspection of the day 807 Gemini im¬ 
ages indicate that the size of the emitting object at the SN location 
is <1."0, or 16 pc at the distance of SN 201 lja. Unfortunately the 
faintness of the object accompanied by poor seeing makes it diffi¬ 
cult to get an accurate FWHM of the object. The average effective 
radii of stellar clusters are ~4 pc, but cluster radii, especially for 


younger clusters, can extend out to 15-20 pc lLarsenll2004t) . mak¬ 
ing it possible that we are indeed seeing a stellar cluster+SN. As 
discussed in more detail below and shown in Figure 8, compari son 
with Starburst99 stellar synthesis models (Lei therer et al,l[l999l) of 
solar metallically give us a decent fit with a 3-6 Myr cluster with the 
same intrinsic extinction of E(B-V) = 1.8. While we cannot defini¬ 
tively rule out the other two options, the late-time luminosity may 
be largely due to the bright parental cluster of SN 201 lja. 

5.4 Progenitor Characteristics 

Mass loss estimates for the progenitor of SN 201 l ja on day 29 and 
113 de rived from X-ray fluxes were presented in lChakraborti et ahl 
d2013l) . Using the narrow Ha emission component present in our 
day 84 spectrum, we estimate the velocity of the surrounding 
medium to be ~180 km s -1 , using the FWHM. This is of course 
the limiting resolution of the spectra, suggesting that this velocity is 
only an upper limit. The true value could be an order of magnitude 
lower and more closely resemble RSG wind speeds. This would 
mean a mass loss of 2 x 10~ 7 - 10~ 6 Mg on day 29 and 3 x 10 -6 
-10~ 5 Mg on day 113, values consistent with BSG and RSG mass- 
loss rates respectively. For comparison HD 168625 is a BSG with 
measured stellar winds of ~ 183 km s -1 and M = 1.2 x 10 -6 Mg, 
dNota et al.ll 1 9961) and SBW1, a near twin of the prog enitor of SN 
1987A , has an estimated mass loss of 3 x 10~ 7 Mf^ ( Smith et al. 
120131) . Therefore as hypothesized by lChakraborti etal.l d20l 3t) . it is 
likely that the increase in X-ray luminosity could be explained by 
the supernova initially expanding into a low-density region then en¬ 
countering a high-density RSG-like environment caused by a RSG 
going through a blue-loop before returning to a RSG phase. If the 
RSG mass-loss had a density enhancement in the equatorial region, 
this could also allow for easy expansion in the polar directions, 
while further increasing the density of the circumstellar ring via 
shocks. 
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There are other possible explanations for the presence of sig¬ 
nificant, asymmetric CSM close to the SN. For example the two 
BSGs listed above, along with Sher25 are all LBV candidates with 
triple ring structures and more importantly, all have a central ring 
surrounding the star. It is possible that the environment created 
around these stars was caused by a violent LBVjtr BSG eruption, 
not density enhancements fro m RSG winds iHendrv et al .1120081 : 
ISmithll2007l : ISrnith et al]|2013l) . It is also possible that the progen¬ 
itor was a more massive RSG that underwent episodic mass loss, 
creating shells a s the faster, less-dense material runs into t he slower, 
denser material ISmith et al.l2009l : lMauerhan et al.l20l 3l) . The ma¬ 
jor factor in all of these scenarios is the inferred initial mass of 
the progenitor. Foran RSG we would expect a progenitor mass < 
25 Mg ( iLevesque et alj|2005t) . whereas for LB Vs we could have 
masses as great as 40-60 Mq. While originally classified as a nor¬ 
mal Type IIP SNe, it is possible that like PTFlliqb (Smith et al. 
2015), SN 201 lja was actually a Type Iln for a very brief stage. 
This could allow for the possibility of a much larger progenitor, 
as it is generally thought that Iln SNe come from LBV or massive 
RSG stars. The explosion date of 2011 December 12 is only an es¬ 
timate, though likely a well constrained one, and it is possible that 
the first observed spectrum could be multiple days older than the 
projected 7 days. PTF1 liqb lost the narrow Iln features ~ 14 days 
after explosion, and any delay in its discovery would have classified 
it as a normal Type II. 

Figure 8 shows the day 807 spectrum fit by both a 3 Myr and 6 
Myr cluster model with solar metallicity produced by Starburst99. 
Without higher S/N and more accurate flux calibration it is impos¬ 
sible to put strict limits on the cluster age and therefore progenitor 
mass. Extensive stellar synthesis modeling and age-dating is be¬ 
yond the scope of this paper, but the spectrum does seem to agree 
with a 3-6 Myr cluster and a progenitor mass > 25 M^. U sing 
high-resolution Keck LRIS spectroscopy, IVinko et~akl (2009) esti¬ 
mate that SN 2004dj had a 10-20 progenitor and res i ded in 
a 10-20 Myr old cluster. More recently Ikuncaravakti et al] d2013l) 
determined a cluster age of 15.6 Myr, corresponding to a mass 
of 14.7 Mg. One notable difference is the lack of Ha emission 
around Sandage-96, whereas pre-explosion Ha imaging of NGC 
4945 from the Danish 1.54 m telescope indicates stro ng emission at 
the location of SN 201 lja llRossa & De ttmar 2003j). Ha emission 
can be a strong indicator of clusterage, and is not commonly found 
in clusters with an age > 8 Myr ^Andrews et al JI20T3I) . Strong Ha 
emission was also seen around SN 2010jl, and archival HST imag¬ 
ing may poin t to a natal cluster < 7 Myr and progenitor mass 
> 30 Mg i j Smith et alj 1201 ll) . It is also interesting to note that 
the BSG Sher25 above may also reside in a ~4 Myr old cluster 
(Melena et apl2008l). _ 

[Chakraborti et al.1 d2013l) conclude the progenitor to SN 
201 lja must be a RSG with a mass > 12 Mq assuming that the 
X-ray variations are only due to a density enhancement by a BSG 
wind and that SN 201 lja was a normal Type IIP SN. The combi¬ 
nation of parent cluster age and Ha retention and 56 Ni mass pre¬ 
sented here, combined with previous X-ray observations all point 
to a larger progenitor mass for SN 201 lja than usually attributed to 
Type IIP SNe. This could be of notable importance since to date, 
observational evidence only suggests RSG progenitor masses of 
Type I IP SNe of up to ~ 16 M^, the so- called “red supergiant prob¬ 
lem” dSmartt et al.ll2009l : ISmartdl2015h . Hypotheses for this mis¬ 
match range from underestimation of dust extinction (Walmswell 
& Eldridge 2012), stellar evolutionary limits (Groh et al. 2013), 
or direct collapse to black holes with little to no optical emission 
(Lovegrove & Woosley 2013, Kochanek 2014). Smith et al. (2011) 


suggest that the high mass RSGs may actually continue to evolve 
into blue or yellow supergiants or Wolf-Rayet stars prior to core 
collapse, forming other CCSNe (IIL. Iln. Ibc). This can be done 
either in a single or binary evolutionary track, and not only allevi¬ 
ates the need to invoke direct collapse to a black hole, but is also 
compatible with the fractional occurrence of non-IIP core-collapse 
events. If in fact the progenitor of the Type IIP SN 201 lja is in the 
20-30 Mq range it may be one of the first definitive examples of 
high-mass RSGs giving rise to a CCSNe. Alternatively, the progen¬ 
itor may be an LBV or a BSG and SN 201 lja may be an example 
of a Type II SNe discovered after a brief Iln phase. 


6 SUMMARY 


We have presented a multi-wavelength analysis of SN 201 lja, an¬ 
other normal Type IIP SN, with not so normal evolution. The short 
optical lightcurve plateau may place SN 201 lja into the realm of 
objects somewhere between IIP and IIL that have been in the spot¬ 
light in recent years. Signatures of dust formation are also seen 
after day 105, when the NIR flux increases, the Ha emission lines 
become attenuated on the red side, and the optical flux decreases. 
At t he same tim e the re is also a noticeable increase in the X-ray 
flux |chakraborti et al. 2013), which is a strong indication of CSM 
interaction. This is further exemplified by the flat-topped hydro¬ 
gen profiles in the optical spectra, particularly the Ha line. The 
dust formation, the flat-topped line profiles, and the increased X- 
ray emission all seem to occur at roughly the same time, which 
paints a picture of the ejecta of SN 201 lja running into dense CSM 
about ~4 months post-explosion. As the ejecta runs into the CSM 
an area is then created between the forward and reverse shock re¬ 
gion where ~6.0 x 10 ~ 5 Mq of dust was quickly formed. It seems 
SN 201 lja is likely another case of dust formation occurring early 
on in a CDS. 

The asymmetry of the line profiles point to a pre-existing CSM 
with a density enhancement inclined 30-45 degrees from edge-on. 
Combined with the possible detection of a young (3-6 Myr) cluster 
at the location of the SN and the large inferred 56 Ni mass of 0.22 
M 0 , it is possible that SN 201 lja had 25-30 Mq LBV or RSG 
progenitor that suffered a massive outburst prior to eruption. Al¬ 
though LB Vs are ge nerally th ought t o be the progenitors of Type 
Iln SNe, recently Smith et all 00151) found that PTFlliqb could 
have been classified as a normal Type IIP SN had it been discov¬ 
ered just days after the initial classification, as the narrow compo¬ 
nent quickly disappeared. This could have been the case with SN 
201 lja as well. If the surrounding CSM was caused by an outburst, 
it would have likely gone undetected due to the extreme extinc¬ 
tion (A„«5.6) in the SN environment. Deep observations with large 
ground- or space-based telescopes as SN 201 lja continues to fade 
will prove essential in putting tighter constraints on the age of the 
cluster and therefore the mass of SN progenitor. 

SN 201 lja seems to share similar properties with the type Iln 
SNe 1998S and 2005ip, the type Ib/c SN 2006jc, and the type II 
SNe 2004dj and 2007od. In all cases interaction with dense CSM 
creates similar observational signatures, and dust formation at early 
times in the CDS regardless of initial SN classification. Particular 
similarities between SNe 201 lja, 2004dj and 2007od indicate there 
is a sub-class of Type IIP SNe that have mass-loss histories capa¬ 
ble of producing pronounced observational signatures normally re¬ 
served for other CCSNe. This only highlights the need for more late 
time observations of even the most seemingly mundane objects, as 
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there seems to be much we do not know about pre-supernova mass 
loss. 
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